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Abstract 

We modify the rules of the self-organized critical forest-fire model in one 
dimension by allowing the fire to jump over holes of < k sites. An analytic 
calculation shows that not only the size distribution of forest clusters but also 
the size distribution of fires is characterized by the same critical exponent 
as in the nearest-neighbor model, i.e. the critical behavior of the model is 
universal. Computer simulations confirm the analytic results. 
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I. INTRODUCTION 



Some years ago, Bak, Tang, and Wiesenfeld introduced the sandpile model which evolves 
into a critical state irrespective of initial conditions and without fine tuning of parameters [|l| . 
Such systems are called self-organized critical (SOC) and exhibit power-law correlations in 
space and time. The concept of SOC has attracted much interest since it might explain the 
origin of fractal structures and 1 / f -noise. Other SOC models e.g. for earthquakes [|2],|3| or the 
evolution of populations |||| have been introduced since then, improving our understanding 
of the mechanisms leading to SOC. Recently, a forest-fire model has been introduced which 
can be viewed as a model for excitable media It becomes self-organized critical when 

time scales are separated f§. In one dimension, the critical exponents could be determined 
analytically, thus proving the possibility of SOC in nonconservative systems 0. 

Analogous to critical phenomena in equilibrium phase transitions, it is expected that the 
values of the critical exponents depend only on few macroscopic properties of the system 
as dimension, conservation laws and symmetries, i.e. the critical behavior of the model is 
universal. Computer simulations of the forest-fire model for different lattice symmetries and 
for a modification with immune trees show indeed universal behavior [|K],[ll[], but so far this 
observation has no analytic foundation. 

In this paper, we show by analytic means that the critical exponents in the one- 
dimensional forest-fire model are universal when the fire is allowed to jump over holes up to 
a given size. In Sec. |l|, we introduce the rules of the model. In Sec. pT| , we give a short 
review of the analytic solution of Ref. |J. In Sec. [IV], we calculate the critical exponent for 
the fire size distribution when the fire is allowed to jump over holes of up to k sites. The 
values of the exponents are confirmed by computer simulations. Finally, we summarize our 
results. 
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II. THE MODEL 



The forest-fire model is a stochastic cellular automaton which is defined on a hypercubic 
lattice with L d sites. In this paper, we consider only the one- dimensional case d — 1. Each 
site is occupied by a tree, a burning tree, or it is empty During one time step, the system 
is parallely updated according to the following rules 

• burning tree — > empty site 

• tree — > burning tree, if at least one neighbor in a distance < k + 1 is burning, 
fc = 0, 1, 2, . . . 

• tree — ► burning tree with probability /, if no neighbor is burning 

• empty site — > tree with probability p. 

Starting with arbitrary initial conditions, the system approaches after a transition period 
a steady state the properties of which depend only on the parameter values. We always 
assume that the lattice is so large that no finite size effects occur. The steady state is 
self-organized critical if the parameters satisfy a double separation of time scales 

f<P<f/p. (1) 

The first inequality means that many trees grow between two lightning strokes and therefore 
large forest clusters and fires occur. The second inequality means that even large forest 
clusters burn down before new trees grow at their edge. Under these conditions, the size 
distributions of forests and fires obey power laws as we shall see below. 

III. ANALYTIC SOLUTION FOR k = 

In the case k — 0, where the fire is stopped by any empty site, i.e. just jumps to 
nearest neighbors, many properties of the model have been derived analytically in 0. Before 
proceeding to general k, we give an illustrative derivation of these results. 
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The mean number of trees destroyed by a lightning stroke is 

s=(f/p)-\l-p)/p, (2) 

where p is the mean forest density in the steady state ||. 

Let n(s) be the mean number of forest clusters of s trees, divided by the number of sites 
L. n(s) will be shown to obey a power law 

n(s) oc s~ T (3) 

for clusters smaller than a cutoff 

Smax OC (f / P )~\ (4) 

besides of possible logarithmic corrections. The probability that lightning strikes a forest 
cluster of size s is proportional to sn(s). Since the fire is stopped by any empty site, the 
size distribution of fires is also proportional to sn(s). 

In order to derive the size distribution of fires, consider a string of n <C p/f sites. This 
string is too short for two trees to grow during the same time step. Lightning does not 
strike this string before all of its trees are grown. Since we are always interested in the 
limit f /p — > 0, the following considerations remain valid even for strings of a very large size. 
Starting with a completely empty state, the string passes through a cycle which is illustrated 
in Fig. |l|. During one time step, a tree grows with probabilty p on any site. After some time, 
the string is completely occupied by trees. Then the forest in the neighborhood of our string 
will also be quite dense. The forest on our string is part of a forest cluster which is much 
larger than n. Eventually that cluster becomes so large that it is struck by lightning with a 
nonvanishing probability. Then the forest cluster burns down, and the string again becomes 
completely empty. (For a rigorous justification of the neglection of lightning strokes on the 
n-string leading to random growth of trees, we refer the reader to Ref. 

This consideration allows us to write down rate equations for the states of the string. 
In the steady state, each configuration of trees is generated as often as it is destroyed. Let 
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P n (m) be the probability that our string is occupied by m trees. Each configuration which 
contains the same number of trees has the same probability. A configuration of m trees is 
destroyed when a tree grows at one of the empty sites, and is generated when a tree grows 
in a state consisting of m — 1 trees. The completely empty state is generated when a dense 
forest burns down. Since all trees on our string burn down simultaneously, this happens 
each time when a given site of our string is on fire. This again happens as often as a new 
tree grows at this given site, i.e. with probability p (1 — p) per time step. We therefore have 
the following equations (which have been derived more formally in |§) 

pnP n (0) =p(l- p), 
p(n — m)P n (m) = p(n — m + l)P n (m — 1) for m/0,n. 

We conclude 

P n (m) — (1 — p)/(n — m) for m < n, 

n—l 

F n (n) = l-(l-/>)£l/(n-m) (5) 

m=0 



i-(i-p)EV 



m. 



m=l 

A forest cluster of size s is a configuration of s neighboring trees with an empty site at 
each end. The size distribution of forest clusters consequently is 

P s+2 (s) l_^p 2 

[) ( s + l)( s + 2)- [ P)S ■ W 



This is a power law with the critical exponent r = 2. The size distribution of fires is 



oc sn(s) oc s 1 . 



There is a characteristic length s max where the power law n(s) oc s~ 2 breaks down. We 
calculate s max from the condition that a string of size n < s max is not struck by lightning 
until all trees are grown. When a string of size n is completely empty at time t — 0, it will 
be occupied by n trees after 

n 

T(n) = (1/p) V] 1/m ~ ln(n)/p 

m=l 



timesteps on an average. The mean number of trees after t timesteps is 

m(t) = n[l — exp(— pt)]. 
The probability that lightning strikes a string of size n before all trees are grown is 

T(n) 

f £ m if) - (f/p)n(\n(n) - 1) ~ (//p)nln(n). 
t=i 

We conclude 

s max ln(s max ) oc p// for large p/ / , (7) 

leading to A = 1. 

Next we determine the relation between the mean forest density p and the parameter 
f/p. The mean forest density is given by 

^ max 

p ~ sn(s) 



vl P) S (« + !)(« + 2) 
1 - p) ln(s max ). 



Thus 



y~r~ p — ln ( s max) ^ ln(p//) for large p/f. (8) 

The forest density approaches the value 1 at the critical point. This is not surprising since 
no infinitely large cluster exists in a one-dimensional system as long as the forest is not 
completely dense. Combining Eqs. (§) and we obtain the final result for the cluster-size 
distribution near the critical point 

,!(S) K (. + !)(.! 2) In.™, fM S < W (9) 
with s max given by Eq. (0). 

The size distribution sn(s) of the fires has also been determined by computer simulations. 
The result is shown in Fig. |^. It agrees perfectly with Eq. (g) in the region s < s max . 



IV. UNIVERSALITY OF THE CRITICAL EXPONENTS 



We now allow the fire to spread to trees up to a distance k + 1 from a burning tree, as 
given by the second rule above. The fire jumps over holes of up to k empty sites, but is 
stopped by holes of more than k sites. Consequently a fire no longer destroys just a single 
forest cluster, but it may also destroy several clusters which are separated by holes of < k 
sites. The size distribution of fires therefore is no longer given by sn(s). In this section, we 
will show that the critical exponent which describes the size distribution of fires is still 1 
indicating its universality. 

We consider a string of n ^> k sites which is empty in the beginning. After t timesteps, 
it contains a hole of size k + 1 with the probability 

n (l - (1 -pY) 2 ~ ne-^ tp 

for small p but large pt. The time after which there are no holes larger than k therefore 
is proportional to ln(n)/p(k + 1). This time becomes very long for large values of n, and 
consequently the forest is very dense at the moment where all holes larger than k have 
disappeared. The critical forest density therefore is still p c = 1. As long as p < 1, there 
is a nonvanishing probability that a hole larger than k occurs, and the fire cannot spread 
indefinitely. 

As in the previous section, we choose n so small that the string is not struck by lightning 
before all empty sites have disappeared, i.e. n < s max . In the limit f /p — > 0, s max diverges, 
and the string can be very large. The dynamics on our string are exactly the same as before 
(Fig. |l]), and Eqs. (§) to (^) are still valid. The size distribution of forest clusters smaller 
than s max remains unchanged, and the critical exponents r = 2 and A = 1 characterizing 
this distribution are universal. Eq. (§) for the forest density also remains the same in the 
limit / /p (besides of a constant which has to be added to the right-hand side of Eq. @ 
but which has already been neglected before, since it is much smaller than ln(p/f)). 

Now we calculate the size distribution F{m) of fires that destroy m trees. For m < s max , 
it can be derived using Eqs. (||). It is proportional to m times the number of configurations 
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which contain to trees with at least k + 1 empty sites at each end and no holes larger than 
k sites between the trees, i.e. 

FyccmV fm+Ww(m) V (m ~ 1)! (io) 

r{m)(xm 2^ ( m+ N+2k+2\ 2^ LI •••L'' 1 ' 

Since we are only interested in the asymtotic power law for large to, this sum can be simpli- 
fied. The main contribution comes from values N <C to since for larger N there are only few 
configurations which contain no holes larger than k (see the consideration at the beginning 
of this section). The probability for a hole of size k + 1 on a string of to + iV sites with 
N «m empty sites is [] 

Am-2) N - k - l /(N -k-l)\ N k+1 
~ m — i)i L Ll L ~ 

~ 1 ; (to - l) N /N\ ~ m k ' 

from which we conclude that the first sum in Eq. flTU| ) has a cutoff for oc m k ^ k+1 \ The 
second sum in Eq. (|ToD counts the number of different configurations of N empty sites on 
m — 1 gaps with the restriction that each gap contains no more than k empty sites. For N 
smaller than the cutoff, the probability for holes larger than k is very small anyway, and we 
are therefore allowed to sum over all configurations of N empty sites on m — 1 gaps which 
gives ( m +!I 2 ~ 2 ). We then obtain 

m k/(k+l) ^ fm+N-2\ m k/(k+i) 2k+1 1 



TO/ TO 



This is a power law with an exponent 1 which is independent of k. We thus have shown that 
not only the exponent for the size distribution of forests but also the exponent for the size 



x This result is valid only for N <C m and is obtained as follows: There are m — l possible positions 
for the hole of size k + 1. The remaining N — k — 1 empty sites can be distributed at random over 
the remaining m — 2 gaps between trees. To obtain a probability, we have to divide by all possible 
configurations of N empty sites on m — 1 gaps. 

2 This is equivalent to the number of configurations of m trees on N + m sites with a tree at both 
ends. 
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distribution of fires is universal. The form of the cutoff functions describing the behavior of 
these distributions on lengths larger than s max , however, is different for different values of k. 

Our computer simulations confirm the analytic result. In Fig. [3] and Fig. [|, the size 
distribution of fires is shown for fire propagation over holes of size 1 and 2. The slope in the 
scaling region is —1 each time. 

V. CONCLUSION 

In this paper we have shown by analytic means that the critical exponents of the SOC 
forest-fire model in one dimension show universal behavior when the range of the interaction 
is changed. This is additionally confirmed by computer simulations. 

In two dimensions, too, computer simulations show that the model is universal under a 
change of the lattice symmetry. It still remains a challenge to prove universality in dimen- 
sions higher than two analytically, e.g. by the renormalization group formalism. 
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FIGURES 

FIG. 1. Dynamics on a string of n = 4 sites. Trees are black, empty sites are white. 

FIG. 2. Size distribution of the fires for f/p = 1/25000 and L = 2 20 . The smooth line is the 
theoretical result, which is valid for cluster sizes < s max - 

FIG. 3. Size distribution of the fires for f/p = 1/8000 and L = 2 20 . The fire is allowed to jump 
over holes of one empty site. 

FIG. 4. Size distribution of the fires for / /p = 1/8000 and L = 2 20 . The fire is allowed to jump 
over holes of two empty sites. 
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